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ABSTRACT 



MS 1603.6+2600 is a high-latitude X-ray binary with a 111 min orbital period, 

,^ . thought to be either an unusual cataclysmic variable or an unusual low-mass X-ray 

binary. In an ASCA observation in 1997 August, we find a burst, whose light curve 

2 ! suggests a Type I (thermonuclear flash) origin. We also find an orbital X-ray modulation 

'^ \ in MS 1603.6+2600, which is likely to be periodic dips, presumably due to azimuthal 

structure in the accretion disk. Both are consistent with this system being a normal 

low-mass X-ray binary harboring a neutron star, but at a great distance. We tentatively 

^ ' suggest that MS 1603.6+2600 is located in the outer halo of the Milky Way, perhaps 

j^ ■ associated with the globular cluster Palomar 14, 11° away from MS 1603.6+2600 on the 

sky at an estimated distance of 73.8 kpc. 

Subject headings: X-rays: stars — stars: individual (MS 1603.6+2600=UW CrB) 



1. Introduction 

MS 1603.6+2600 was first discovered in the course of the Extended Medium Sensitivity Survey 
(Gioia et al 1990). With only 51 detected photons in a 2112 s exposure, this detection was useful 
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only in providing the position and a flux estimate of ~ 10~^^ ergscni^^s^^ (0.3-3.5 keV). Morris 
et al (1990) identified the optical counterpart, subsequently designated UW Coronae Borealis, a V 
~19.7th mag object. It was shown to be an eclipsing binary: measurements of 10 eclipse timings 
over a 3-day period has enabled Morris et al (1990) to derive a 111.0 min orbital period. This 
immediately establishes MS 1603.6+2600 as a low-mass, compact binary. To fit in such a short 
period binary, the secondary must be a low mass star; and, to be a significant X-ray source, the 
primary must be an accreting compact object — a white dwarf, which would make the system a 
cataclysmic variable (CV), a neutron star, or perhaps a black hole, either of which would make it 
a low mass X-ray binary (LMXB). Given the high galactic latitude of this object (6^^ ~ 47°), it 
has to be at a large distance (d ^ 10 kpc) to be a normal, bright LMXB. On the other hand, the 
optical spectrum and the high X-ray to optical flux ratio (fx/fopt -^15) are unlike those of normal 
CVs. 

Further optical observations (Vilhu et al 1993; Hakala et al 1998) have revealed highly variable 
orbital light curves. Report of a 112.5 min period, obtained by Vilhu et al (1993) from Fourier 
analysis, would make it a system with multiple periods, if conflrmed. A ROSAT PSPC observation 
performed on 1991 August 26 (Hakala et al 1998) did not resolve the nature of MS 1603.6+2600. 
Hakala et al (1998) discussed the possibility that MS 1603.6+2600 may be an LMXB with an 
accretion disk corona (ADC), in which we only observe a small fraction of X-rays scattered in the 
ADC surrounding the inner disk and the neutron star. They also considered the possibility that it 
might be a transient LMXB in quiescence. These scenarios are motivated, in part, by the desire to 
keep MS 1603.6+2600 within a few kiloparsecs of the Galactic plane. 

We have obtained an ASCA observation of MS 1603.6+2600 in an effort to clarify the nature 
of this object (see also our preliminary report in Mukai et al (1999)). We describe the observation 
in §2, the results in §3, and discuss the implications in §4. 



2. Observation 

MS 1603.6+2600 was observed with ASCA (Tanaka et al 1994) from 1997 August 30 at 
20:25 UT to 1997 August 31 at 13:52 UT. Of the 4 instruments on-board ASCA, the 2 GIS's were 
operated in the standard PH mode. For the CIS data, we have selected intervals when the satellite 
was outside the South Atlantic Anomaly (SAA), the attitude control was stable and within 0.02° 
of the target, and the line-of-sight was > 5° above the Earth limb. We have also applied a standard 
selection expression combining monitor count rates and geomagnetic cut-off rigidity (COR) to 
exclude time intervals of high particle background, obtaining ~25 ksec of good on-source data. We 
extracted the source spectra and light curves using 6 arcmin radius circular extraction regions, and 
used ~13 arcmin radius regions around the center of the detectors, excluding regions within ~7.5 
arcmin of the source, as background (no sources were detected in the background regions both for 
CIS and for SIS). The net count rates were ~0.06 ctss~^ in GIS-2 and ~0.07 ctss~^ in GIS-3. 
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Fig. 1. — (Upper panel) ASCA light curve of the burst, in 16 s bins. We show the total count rate 
for 2 SIS (0.48-10 keV) and 2 GIS (0.7-10 keV) instruments. (Lower panel) The background light 
curve during the same interval. 



For the 2 SIS instruments, we have applied the same SAA, attitude, and elevation constraints, 
and further excluded data taken within 10° of the bright Earth limb. We have selected low back- 
ground time intervals using COR and PIXEL monitor count criteria, and also excluded data taken 
within 32 seconds of day/night and SAA transitions, when the on-board dark frame calculations 
are suspect. These resulted in ~25 ksec of useful SIS data. We used an extraction region of --^4.5 
arcmin radius for the SIS-0, using the entire chip, excluding regions within ~5.5 arcmin of the 
source, as background, with a net rate of '^0.12 ctss^-*^. For SIS-1, we were forced to use a smaller 
('^3.5 arcmin) source region to stay entirely within the active chip, obtaining a net rate of ~0.08 
ctss^^. 

For spectral analysis, we combined data from like instruments to form a GIS spectrum and an 
SIS spectrum. The latter is affected by the secular change in SIS performance at low energies; this 
can be phenomenologically characterized as a spurious excess absorption of order Nh = 1.0 x 10^^ 
cm~^. For the light curve analysis, we combined data from all 4 instruments, both using the entire 
ASCA passband, and in two sub-bands, one below 2 keV (0.48-2.0 for SIS and 0.7-2.0 keV for GIS) 
and the other above 2 keV (up to 10 keV). 



3. Results 



3.1. Detection of a burst 

We present a 16-s bin light curve of a striking, flare- like event in Figure 1. This event, detected 
at 21:45:50 UT on 1997 August 30, is seen in all 4 individual source light curves, and in none of 
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Fig. 2.— The ASCA light curves of MS 1603.6+2600 above and below 2 keV, folded on the orbital 
ephemeris of Morris et al (1990), in 32 bins per cycle, plotted twice for clarity. A hardness ratio 
plot is shown on the top panel in 16 bins per cycle. 



the 4 background light curves. We are therefore confident that this is a real event, presumably due 
to MS 1603.6+2600. The rise is rapid, and the decay lasts for about 1 minute before blending into 
the (noisy) persistent emission. The peak count rate (~ 2.6 ctss~^ in 4 detectors) is 8 times the 
persistent emission. 

We have attempted a more detailed analysis of this event, but this has not been fruitful because 
there are only ~60 net counts from this event. In the remainder of this paper, we will refer to this 
event as a "burst." The resemblance with a Type I X-ray burst (i.e., a thermonuclear runaway on 
the surface of an accreting neutron star; see Lewin et al (1995) for a review) is suggestive, but not 
conclusive. 



3.2. Orbital modulation 



Next, we removed the burst from our light curves and performed Fourier and folding analyses. 
We find a strong periodicity at 116+6 min, the large uncertainty being due to the limited duration 
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Fig. 3.— The archival ROSAT FSPC Ught curve of MS 1603.6+2600, obtained on 1991 August 26. 
The data have been divided into two sections, and folded on the orbital period with 20 bins. 



of the observation and data gaps. We interpret this signal as reflecting the orbital modulation, and 
proceed below by folding the data on the orbital ephemeris of Morris et al (1990). Note, however, 
that this ephemeris does not define phase 0.0 accurately. No other periodicities are apparent in the 
data in the 200-20,000 s range. 

We present the folded orbital light curves, in high (>2 keV) and low (<2 keV) energy bands, 
as well as the hardness ratio, in Figure 2. A clear orbital minimum, lasting ~15% of the orbit and 
reaching ~60% of the flux level outside the minimum, is seen in the low energy light curve. It is less 
clear in the high energy curve, which is reflected in the slight increase in the hardness ratio at these 
phase intervals. Inspection of individual orbital cycles shows that the minimum is always present 
during the same phase in broad outline, though their details vary. However, given the statistics, 
we cannot rule out statistical fluctuation as the cause of the cycle-to-cycle variations. Outside the 
minimum, the light curves are flat, with apparently random fluctuations superimposed. 

This behavior is rather different from that reported by Hakala et al (1998) in their 1991 
ROSAT PSPC observation. To investigate this further, we have retrieved the ROSAT data from 
the HEASARC and re-extracted the light curves. At first glance, we merely confirmed the results 
of Hakala et al (1998). However, after dividing the ROSAT data into two halves of roughly 
equal durations (~12 hrs), we find a pattern of orbital minimum for the first half of the ROSAT 
observation (Figure 3). This minimum is absent during the second half; instead, the folded light 
curve shows a flaring-like behavior. 
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Fig. 4. — ASCA GIS spectra of MS 1603.6+2600 during orbital maximum (pluses) and minimum 
(diamonds). The top panel shows the data and the model (see text for details), while the bottom 
panel shows the data to model ratio. 



3.3. Spectral analysis 



Despite the improvement over the ROSAT PSPC spectrum, we have not been able to draw 
a definite conclusion from the average ASCA spectrum of MS 1603.6+2600. It can be fit with a 
power law, a bremsstrahlung, or a Comptonization model: discrimination among these continuum 
models proved impossible because high signal-to-noise ratio was not achieved over a sufficiently wide 
bandpass. No discrete features are obvious. Regardless of the continuum model, the 0.7-10 keV 
flux of MS 1603.6+2600 during the ASCA observation was ~4.0 x 10"^^ grgscm^^ g-i (^43 ^ 10^^ 
[d/10 kpc]^ ergss~^). Ignoring the difference in spectral shapes, the --^8 fold increase in count rate 
at the peak of the burst implies a peak luminosity of '^3.8 x 10"^^ [d/10 kpc]^ ergss"^). 



The average ASCA flux value indicates a significant long-term variability in MS 1603.6+2600. 
The Einstein IPC 0.3-3.5 keV flux of 1.14 xlO~^^ ergscm^^s^^ (Morris et al 1990) is roughly 
a factor of 4 lower than the ASCA value, although both the bandpass and the assumed spectral 
shape are different. Back prediction of Einstein count rate based on ASCA flux and the best-fit 
bremsstrahlung model also points to a factor of ~4 change in flux. Similarly, MS 1603.6+2600 
appears to have been fainter by a factor of ~7 during the i? 05^ T observation, compared to 1997 
August. 

We also compared the spectrum during the orbital minimum to that of the maximum, by 
simultaneously fitting the two using the same bremsstrahlung model (identical kT and normaliza- 
tion). For the maximum, we assumed no absorption, since the interstellar column, of the order 
of 10^*^ cm^^ according to the ROSAT PSPC spectral fit (Hakala et al 1998), is too small to 
be detectable with ASCA. For the minimum spectrum, we have added an absorption component. 
Either a simple absorber of Nh = 4.6 x 10^^ cm^^ (x^=70.6 for 140 PHA bins and 4 degrees of 



freedom) or a partial covering absorber of Nh = 1.3 x 10 cm and a covering fraction of 37% 



'0 



(x^=64.8 for 140 PHA bins and 5 degrees of freedom) can be used to fit both spectra. The latter 
is shown for the GIS data in Figure 4. 



4. Discussion: the Nature and the Location of MS 1603.6+2600 

We have observed an orbital minimum in the ASCA light curve which is broad and partial. 
Our re-analysis of the ROSAT PSPC light curve shows that a similar minimum was present during 
the first half of the 1991 observation, but disappeared during the second half. 

In the ADC model that we favored in our earlier report (Mukai et al 1999), this minimum 
is caused by a partial eclipse of an extended X-ray source, whose size is a significant fraction of 
the binary separation. An ADC source is an LMXB seen at a high inclination angle. The central 
X-ray emitting region is blocked from our direct view by the accretion disk. A fraction (1-10%) is 
scattered into our line of site by the ADC above the accretion disk. Although this interpretation 
has several attractive properties, MS 1603.6+2600 does not show the smooth, quasi-sinusoidal 
modulation which is seen in the prototype ADC source, X1822— 372 (Hellier & Mason 1989). The 
orbital light curves of ADC sources generally appear to be stable from one cycle to the next, while 
the ROSAT data show MS 1603.6+2600 to be otherwise. In addition, the spectrum of X1822-372 
is softer during the partial eclipse, which is the reverse of what we observe in MS 1603.6+2600. 
Therefore, the behavior of MS 1603.6+2600 differs in detail from those of well-established ADC 
sources such as X1822— 372. Moreover, the ingress into the orbital minimum in MS 1603.6+2600 
(Figure 2) takes place in less than one phase bin (l/32th of the orbital cycle). It is difficult to 
construct a geometry of an extended X-ray source that leads to a broad, partial eclipse yet with 
such a relatively rapid ingress. 

On the other hand, the behavior of MS 1603.6+2600 is similar to that of dipping LMXBs such 
as XB 1916—053 (Homer et al 2001). In dippers, we believe that the central X-ray source (the 
accretion disk, the neutron star, or both) is occulted by azimuthal structure of the accretion disk. 
The preferential location of the structure in binary frame defines the envelope of phases during 
which dips occur. The dip durations and depth change from epoch to epoch, sometimes within a 
single observation. Spectrally, dips are deeper at low energies, although simple absorption models 
often do not fit the data, instead requiring a partial covering absorber. In all these respects, the 
ASCA data of MS 1603.6+2600 are consistent with the dipper interpretation. We therefore favor 
this interpretation over the ADC model. 

As a dipper, MS 1603.6+2600 must possess an accretion disk. This immediately excludes the 
AM Her-type magnetic CVs (in which the magnetic field is so strong as to prevent the formation of a 
disk) from consideration. It is still possible that MS 1603.6+2600 is a dipping CV with an accretion 
disk. For example, U Gem, the prototype dwarf nova, is a dipper (Szkody et al 1996). However, 
all the objections raised to date remain valid, most notably that the fx/fopt ratio is too high: U 



Gem is ~1,000 times brighter than MS 1603.6+2600 in the optical while only ~4 times brighter in 
the ASCA band. Magnetic CVs of the intermediate polar type, which has partial accretion disk, 
also show X-ray dips (Hellier et al 1993). However, their defining characteristic is a strong spin 
modulation in the X-rays, which we do not see. Thus, if MS 1603.6+2600 is a CV, it must be an 
unusual example. 

Moreover, a burst-like event like the one we have seen in MS 1603.6+2600 has never been 
reported in any X-ray observations of CVs. We therefore assume, as a working hypothesis, that 
MS 1603.6+2600 is an LMXB containing a neutron star and the burst is a Type I event. This 
suggest MS 1603.6+2600 is at a great distance, since the typical persistent luminosity of a dipper is 
a few times 10^^ ergss"^ and the typical burst peak luminosity is a few times 10^^ ergss~^. Given 
the observed flux values, this would require d^ 10 kpc. 

Is this a serious weakness of this model? Is there anything in the Galactic halo out to, say, 
d ~50 kpc? To answer these questions, we have conducted a simple search for globular clusters 
near the direction of MS 1603.6+2600. The nearest match was found to be Palomar 14 (Sarajedini 
1997), --^11° away on the sky at an estimated heliocentric distance of 73.8 kpc. 

This suggests possible, though speculative, scenarios as to how MS 1603.6+2600 might have 
formed in the outer halo. First, it may have formed in Palomar 14 itself and escaped. Of the -^150 
Galactic globular clusters, 12 harbor a luminous X-ray binary (Deutsch et al 2000). The same 
stellar encounters that form these LMXBs eject a significant number of resultant binaries. One 
calculation (Portegies Zwart et al 1997) suggest the fraction of the escaped binaries may be about 
a third of those that have survived in the globular cluster, thus we expect several "non-globular- 
cluster" LMXBs to have escaped from the parent cluster. Similar escapees from a cluster much 
closer to the Galactic center would have been assimilated into the Galactic bulge, and its origin 
would have become obscured. If MS 1603.6+2600 was ejected from Palomar 14 at a velocity of 10 
kms~^, it would have taken it -^1 Gyr to move 10 kpc, which is the minimum distance implied by 
the current angular separation of 11°. Thus it appears possible to have a reasonable combination 
of space velocity and lifetime of MS 1603.6+2600 to explain its current location, for an assumed 
origin in Palomar 14. 

Another speculative scenario is related to the possible origin of Palomar 14, and many other 
outer halo globular clusters, in dwarf elliptical satellites of the Milky Way (van den Bergh 2000). 
Several outer halo clusters are associated with the Sagittarius dwarf. Similarly, Palomar 14 may be 
associated with an undiscovered extant satellite, or one that has been tidally disrupted. In either 
case, MS 1603.6+2600 could be a member of this hypothetical satellite galaxy. 

If we assume a physical association between MS 1603.6+2600 and Palomar 14, at a distance of 
d rsj 70 kpc, then the average luminosity during ASCA observation was --^2.4 x 10^^ ergss"^ and the 
peak burst luminosity was ~1.9 x 10^^ ergss"^, both of which are normal for a bursting, dipping 
LMXB, although still with an unusually low ix/^opt ratio. Evolutionary scenario la of Ergma & 
Vilhu (1993), similar to short period CVs, predicts a similar luminosity and bursting behavior, and 



hence fits our outer halo interpretation. However, this may be a coincidence if MS 1603.6+2600 is 
indeed formed in Palomar 14, since globular cluster LMXBs form in a distinct manner from field 
LMXBs. 

If our interpretation is correct, MS 1603.6+2600 is a close analog of XB 1916—053, including 
the strong variabilities of the dip morphology. In fact, the second optical period obtained via 
Fourier transform (Vilhu et al 1993) may be real, indicating that MS 1603.6+2600 is a permanent 
superhumper (Haswell et al 2001), another characteristic of XB 1916—053. If so, it would be 
possible to estimate the mass ratio of MS 1603.6+2600 using the fractional period excess. 



5. Summary and Future Prospects 

We have presented our ASCA observation of the unusual high-latitude X-ray binary, MS 1603.6+2600. 
We have detected a burst and an orbital modulation consistent with a dipping behavior, and pre- 
sented an outer halo LMXB model as perhaps the least contrived interpretation. 

Future X-ray observations can test our interpretation. In particular, the detection of other 
bursts with better sensitivity is necessary to confirm the Type I interpretation spectroscopically. 
Moreover, if we are lucky enough to catch a radius expansion burst, a reliable distance and other 
parameters can be derived. The ingress and egress to the orbital minimum should be investigated 
further, so that we can can conclusively choose between the dipper and the ADC interpretations. 
At the same time, further optical observations should clear up the issue of possible multiple peri- 
odicities. 

Pending such observations, we tentatively conclude that MS 1603.6+2600 may be the most 
distant Galactic X-ray binary known. Moreover, this may be an escaped globular cluster LMXB, 
something we cannot confidently observe near the Galactic center, or may be a relic of the hierar- 
chical formation of the Milky Way halo. 

This research has made use of data obtained from the High Energy Astrophysics Science 
Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. RW is 
supported by NASA through Chandra Postdoctoral Fellowship grant PF9-10010 awarded by CXC, 
which is operated by SAS for NASA under contract NAS8-39073. 
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